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Abstract. The model to simulate jet quenching effect in ultrarelativistic heavy ion 
collisions is presented. The model is the fast Monte-Carlo tool implemented to modify 
a standard PYTHIA jet event. The model has been generalized to the case of the 
"full" heavy ion event (the superposition of soft, hydro-type state and hard multi-jets) 
using a simple and fast simulation procedure for soft particle production. The model 
is capable of reproducing main features of the jet quenching pattern at RHIC and is 
applyed to analyze novel jet quenching features at LHC. 



PACS numbers: 25.75.-q, 12.38.Mh, 24.85. +p 
1. Introduction: 

One of the important tools for studying the properties of quark-gluon plasma (QGP) 
in ultrarelativistic heavy ion collisions is the analysis of a QCD jet production. 
The medium-induced energy loss of energetic partons, "jet quenching", should be 
very different in the cold nuclear matter and QGP, resulting in many observable 
phenomena p]. Recent RHIC data on high-p^ particle production at ^/s = 200 A 
GeV are in agreement with the jet quenching hypothesis [2]. At LHC, a new regime 
of heavy ion physics will be reached at ^/snn — 5.5A TeV where hard and semi-hard 
particle production can stand out against the underlying soft events. The initial gluon 
densities in Pb+Pb reactions at LHC are expected to be much higher than those at 
RHIC, implying a stronger partonic energy loss, observable in new channels. 

In the most of available Monte-Carlo heavy ion event generators the medium- 
induced partonic rescattering and energy loss are either ignored or implemented 
insufficiently. Thus, in order to analyze RHIC data on high-p^ hadron production 
and to test the sensitivity of LHC observables to the QGP formation, the development 
of adequate and fast Monte- Carlo tool to simulate the jet quenching is necessary. 



2. Physics model and simulation procedure 

The detailed description of physics model can be found in our recent paper [3]. The 
approach bases on an accumulating energy loss, the gluon radiation being associated 
with each parton scattering in the longitudinally expanding quark-gluon fluid and 
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includes the interference effect using the modified radiation spectrum dE/dl as a function 
of decreasing temperature T. The basic kinetic integral equation for the energy loss AE 
as a function of initial energy E and path length L has the form 

v ' 

where / is the current transverse coordinate of a parton, dP/dl is the scattering 
probability density, dE/dl is the energy loss per unit length, A is in-medium mean 
free path. The collisional loss in high-momentum transfer limit, radiative loss in 
BDMS approximation [I], simple Gaussian parameterization of gluon emission angle 
distribution and realistic nuclear geometry are used. The QGP proper formation time 
To and initial temperature To are model parameters. 

The model has been constructed as the Monte-Carlo event generator PYQUEN 
(PYthia QUENched) and is available via Internet [5]. The event-by-event simulation 
procedure includes the generation of the initial parton spectra with PYTHIA [6] and 
production vertexes at given impact parameter, rescattering-by-rescattering simulation 
of the parton path length in a dense zone, radiative and collisional energy loss per 
rescattering, final hadronization with the Lund string model for hard partons and in- 
medium emitted gluons. 

The full heavy ion event is simulated as a superposition of soft hydro-type state 
and hard multi-jets. The simple approximation [3] of hadronic liquid at "freeze-out" 
stage has been used to treat soft part of the event. Then the hard part of the event 
includes PYQUEN multi-jets generated according to the binomial distribution. The 
mean number of jets produced in AA events at a given impact parameter is a product 
of the number of binary NN sub-collisions and the integral cross section of hard process 
in pp collisions with the minimal transverse momentum transfer p™ in . The extended in 
such a way model has been also constructed as the fast Monte-Carlo event generator [TJ. 
Note that ideologically similar approximation has been developed in [8]. 



3. Jet quenching at RHIC 

In order to demonstrate the efficiency of the model, the jet quenching pattern in Au+Au 
collisions at RHIC was considered. Model parameters has been selected from the fit of 
rj— and pt— spectra at RHIC for different event centralities (see [3j for details). In 
particular, the nuclear modification of the hardest domain of p^-spectrum was used to 
extract initial QGP conditions: To = 500 MeV and tq = 0.4 fm/c. Figure [U shows that 
our model reproduce pt— and centrality dependences of nuclear modification factor 
Raa (determined as a ratio of particle yields in AA and pp collisions normalized on 
the number of binary NN sub-collisions) for 7To's measured by PHENIX [9] quite well. 
Another important tool to verify jet quenching is two-particle azimuthal correlation 
function C(A<p) - the distribution over an azimuthal angle of high-p^ hadrons in the 
event with 2 GeV/c < px < p^ relative to that for the hardest "trigger" particle with 
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p£ lg > 4 GeV/c. Figure [2] presents C[A(p) in pp and in central Au+Au collisions (data 
from STAR [10J). Clear peaks in pp collisions at A(p = and Atp = it indicate a typical 
dijet event topology. However, for central Au+Au collisions the peak near 7r disappears. 
It can be interpreted as the observation of monojet events due to the absorption of one 
of the jets in a dense medium. Figure [2] demonstrates that measured suppression of 
azimuthal back-to-back correlations is well reproduced by our model. 
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Figure 1. The nuclear modifi- 
cation factor Raa for neutral pi- 
ons in Au+Au collisions for two 
centrality sets. The points are 
PHENIX data, histograms are 
the model calculations. 



Figure 2. The azimuthal two- 
particle correlation function for 
pp and for central Au+Au colli- 
sions. The points are STAR data, 
dashed and solid histograms are 
the model calculations for pp and 
Au+Au events respectively. 



We leave beyond the scope of this paper the analysis of such important RHIC 
observables as the azimuthal anisotropy and particle ratios, which are sensitive to the 
soft physics. In order to study them, a more careful treatment of low-pr particle 
production than our simple approach is needed [TT] (the detailed description of space- 
time structure of freeze-out region, resonance decays, etc.). 



4. Jet quenching at LHC 



The developed model was applied to analyze various novel features of jet quenching at 
the LHC (see [T5J ESI EH] for details). Let us just to enumerate some main issues. 

Jets tagged by leading hadrons. The relation between in-medium softening jet 
fragmentation function (measured with leading hadrons) and suppression of jet rates due 
to energy loss out of jet cone, was analyzed in the work [T2]. The specific anti-correlation 
between two effects allows one to probe parton energy loss mechanism (small-angular 
radiation versus wide-angular radiation and collisional loss). 
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Jets induced by heavy quarks. The possibility to observe B-jets tagged by an 
energetic muon in heavy ion collisions was analyzed in the work |13j . The significant 
softening b-jet fragmentation function due to b-quark energy loss is predicted. 

Z/^+jet production. The channel with dimuon tagged jet production in heavy 
ion collisions was analyzed in the work [14] . We have found that the medium-induced 
partonic energy loss can result in significant P^r-imbalance between fi + n~ pair and a 
leading particle in a jet, which is quite visible even for moderate loss. 



5. Conclusions 



The method to simulate jet quenching in heavy ion collisions has been developed. 
The model is the fast Monte-Carlo tool implemented to modify a standard PYTHIA 
jet event. The full heavy ion event is obtained as a superposition of a soft hydro- 
type state and hard multi-jets. The model is capable of reproducing main features of 
the jet quenching pattern at RHIC (the px~ dependence of the nuclear modification 
factor and the suppression of azimuthal back-to-back correlations). The model was also 
applied to probe jet quenching in various new channels at LHC energy: jets tagged by 
leading particles, b-jets, dilepton-jet correlations. The further development of the model 
focusing on a more detailed description of low-pr particle production is in the progress. 
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